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Three low-valence layered nickelates with general formula LasN^Or-^ (£=0.0, 0.5, 1) are studied 
by ab initio techniques. Both the insulating and metallic limits are analyzed, together with the 
compound at the Mott transition (5= 0.5; Ni 2+ ), that shows insulating behavior, with all Ni atoms 
in a S=l high-spin state. The compound in the 6= 1 limit (LasN^Oe), with mean formal valence 
Ni 1-5+ and hence nominally metallic, nevertheless shows a correlated molecular insulating state, 
produced by the quantum confinement of the Ni02 bilayers and the presence of mainly d z 2 bands 
(bonding- antibonding split) around the gap. The metallic compound shows a larger bandwidth of 
the e g states that can sustain the experimentally observed paramagnetic metallic properties. The 
evolution of the in-plane antiferromagnetic coupling with the oxygen content is discussed, and also 
the similarities of this series of compounds with the layered superconducting cuprates. 



I. INTRODUCTION 

Three types of metal-insulator transitions have 
been observed in nickel oxides. The first one oc- 
curs in the RMO3 perovskites, by changing the R 
cation. LaNiOs is a paramagnetic metal with en- 
hanced correlations, but as the size of the lanthanide 
cation is increased, the series of RMO3 (with R=Y 
and from Sm to Lu in the lanthanide series) are in- 
sulators that undergo a temperature induced metal- 
insulator transition, 1 with the metal-insulator tran- 
sition temperature (accompanied by a transition to 
a low-temperature antiferromagnetic state) grows as 
the R-0 bond length increases. 2 This transition is 
accompanied of a charge-ordering (or at least charge 
disproportionation) phenomenon, that has been ob- 
served in various members on the insulating side of 
the transition. 3 

A second type of metal-insulator transition has 
been observed very recently in ultrathin films of 
LaNiOs grown as multilayers on LaAlO^ and also as 
thin films £ In the multilayered case, spectroscopic 
studies showing a double peak structure developing 
in the insulating cases suggest the appearance of 
charge-order in the Ni layer, similar to what hap- 
pens in bulk RMO3 perovskites for certain types of 
rare earth other than L a, as discussed above. It is 
likely that charge order is the origin of the insulating 
behavior in the case of thin films as well, but no spec- 
troscopy has been performed yet to shed light on the 
cause of the insulating behavior observed. A similar 
phenomenon of charge disproportionation has been 
reported very recently 6 at the interface between two 
SrTiOa thin films when a monolayer of RO is grown 
at the interface, where R can be various lanthanides. 
It turns out that for the La case the interface is 
metallic, but for Y and the other lanthanides it is 
insulating, due to a similar phenomenon of charge 
ordering as in perovskite nickelates. 

The third type of metal-insulator transition that 



has been experimentally observed in nickelates ap- 
pears in the layered Ruddelsden-Popper series for 
n=2 and 3, i.e. La3Ni208-5 and I^NisOio-^. 7 
La4Ni30io and La3Ni20y are both paramagnetic 
metals* 7 - - — However, at a certain value of J, 
a doping-induced metal-insulator transition takes 
place. Experimentally, various doping levels have 
been reached, up to La3Ni206.35 and La4Ni30s.76, 
showing this metal-insulator transition. 7 In princi- 
ple, a Mott-transition could take place at the Ni 2+ 
valence (La3Ni2 06.5 or L^NisOg), making those 
systems already insulating. This has to do with the 
formation of a Ni 2+ :d 8 electronic state, which in an 
octahedral environment would form a S=l state with 
an e^ 2 electronic configuration. The e p t band will be 
fully occupied leading to a gap opening if the Hund's 
coupling is larger than the bandwidth of the e g bands 
and no big distortion from octahedral environment 
occurs that would split the otherwise degenerate e^ 
bands. 

The situation would be reminiscent of what hap- 
pens on La2Ni04, which also, by electron count, has 
nominally Ni 2+ cations. The material shows insulat- 
ing behavior at low temperature, and also a peculiar 
type of metal-insulator transition at high tempera- 
tures, that has been explained by Goodenough 10 as 
a splitting in the semi- itinerant d x 2_ y 2 band that 
opens up a gap, with the d z 2 behaving always as lo- 
calized electrons, according to his picture. How far 
the insulating region reaches in 5 values will depend 
on how quick metallicity is reached when the insu- 
lating system is doped with either holes or electrons. 
Experiments show that even La3Ni2 06.84 (formal va- 
lence Ni 2,34+ ) is an insulating compound^ 7 - suggest- 
ing that the insulating phase is difficult to destroy by 
doping, more difficult than in the case of cuprates, 
probably due to the relatively small polaron size in 
doped nickelates. Below, we will discuss the dis- 
tinctions compared to layered cuprates that might 
help explaining this difference in behavior. The use 
of these low- valence layered nickelates has been re- 



cently suggested as a possible way to mimic the elec- 
tronic structure of superconducting cuprates j 11 ^ 2 

Other Ni 2+ insulating compounds have been 
found, e.g. LaNi02.5, with a superstructure based 
on the perovskite unit cell^i 3 - When lowering the va- 
lence below Ni 2+ , the layered nickelates that have 
been synthesized remain insulating. Experiments on 
recently synthesized La4Ni308— ~— show it is an in- 
sulator with antiferromagnetic (AF) order that leads 
to a peculiarity in the susceptibility at high fields. 
In the case of I^NiaOg, an unusual type of molec- 
ular insulating state has been used to understand 
its electronic and magnetic properties^ where the 
quantum confinement of the MO2 trilayers leads to a 
bonding-antibonding splitting that creates an insu- 
lating state with a correlated molecular origin. Be- 
low, we will see how a similar picture applies to 
LasNi206, which is also insulating, but no AF or- 
der has been observed down to 4 K*2°- 

Even lower Ni valence layered compounds exist 
and again a crossover from insulating to metallic 
state occurs as valence is reduced. The infinite-layer 
LaNi02 has been synthesized, with single-valent 
Ni + :d 9 cations, 21 showing a paramagnetic metal- 
lic behavior. Electronic structure calculations show 
the origin of the metallicity of this compound and 
study the differences with the isoelectronic parent 
compounds of the high-temperature superconduct- 
ing cuprates. 22 



II. STRUCTURE 

We have studied the n= 2 members of the 
Ruddelsden-Popper 23 series La n +iNi n 03 n +i, with 
reduced oxygen content, to explore the regime of Ni 
valencies around 2+. These structures are charac- 
terized by the presence of two MO2 layers separated 
by a block of fluorite La/02 /La layer that produces 
the confinement of the MO2 bilayer in the structure 
(see Fig. [1]) and determines the valence. 

As the oxygen content is varied, so does the Ni co- 
ordination. In the insulating phase La3Ni2 06 (right 
panel of Fig. [1]), the Ni environments are square pla- 
nar. However, in the metallic La3Ni2 07, every Ni is 
in an octahedral oxygen environment (left panel of 
Fig. [1]). On going from one oxygen content in the 
metallic limit to the other in the insulating side of 
the transition, apical oxygens need to be removed, 
that would otherwise be in the Ni bi-layer. In the in- 
termediate case, La3Ni2 06.5, we see that the reduc- 
tion of oxygens leads to a mixture of octahedral and 
square-pyramidal environments. Structures used for 
the calculations in the insulating and metallic limits 
were taken from Refs. I24U25I . For the relaxation of 
the intermediate structure, we have considered the 



possible apical oxygen removals and compared their 
total energies, leading to the structure presented in 
the central panel of Fig. [TJ In our calculations, we 
have studied an ordered alternation of octahedra and 
square pyramids. Such an ordering need not happen 
in the real compound, but it is a specific way to 
calculate its electronic structure using a simple unit 
cell. One should keep in mind that the positioning 
of two neighboring square pyramids along the c-axis 
would change the d 2 2-d 2 2 a-bonding that has cru- 
cial consequences in the electronic structure of the 
compound (see below). When analyzing the possible 
bonding-antibonding splittings produced by such a 
structure, it is important to keep in mind that a ran- 
dom distribution of square pyramids and octahedra 
would modify them. 



III. COMPUTATIONAL DETAILS 

Our electronic structure calculations were per- 
formed within density functional theor y 26 ! 27 using 
the all-electron, full potential code wien2k 28 based 
on the augmented plane wave plus local orbital 
(APW+lo) basis set. 2 - The generalized gradient 
approximation 3 - - (GGA) was used for the structure 
optimization of the 5= 0.5 compound for various 
ways to obtain that oxygen content starting from 5= 
0. To deal with strong correlation effects we apply 
the LDA+U scheme 31,32 including an on-site repul- 
sion U and Hund's coupling J for the Ni 3d states. 
Results that we present are not dependent on the 
specific values of U and J within a reasonable range, 
and we report results with U = 4.75 eV, J = 0.68 
eV, values very similar to those determined from 
constrained density functional calculations 16 for a 
similar layered nickelate. 



IV. LA 3 NI 2 6 

La3Ni2 06 has been studied experimentally re- 
cently, and also using ab initio calculations. 20 The 
system shows largely insulating behavior that can 
be described by a variable range hopping model. No 
magnetic ordering is found down to 4 K. The calcu- 
lations presented in that paper report only metallic 
results, hence they do not provide an understand- 
ing of the observed insulating behavior of the com- 
pound. In our recent paper, 19 we describe the origin 
of the insulating behavior of the similar compound 
La4Ni308 as being due to the gap opening that is 
produced by a bonding-antibonding splitting of the 
d z 2 bands closest to the Fermi level. This comes 
about due to the spatial confinement of the MO2 
trilayers in the compound. In the case of La3Ni206, 
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FIG. 1: (color online) Structure of the layered nickelates, removing oxygen from left to right: LasNi207, the metallic 
limit on the left side; La4Nis06.5, the Mott-transition oxygen content in the middle; and La4Nis06, the molecular 
correlated insulating limit on the right. Large (green) sphers denote La, oxygen lies at the vertices of the octahedra, 
pyramids, and squares. 



the structure presents a MO2 bilayer, but the situ- 
ation becomes similar to La4Ni30s. Figure [2] shows 
the band structure of the compound, with the dif- 
ferent Ni e g states highlighted, calculated within the 
LDA+U scheme. Due to the square planar environ- 
ment of the Ni 1,5+ cations (assuming a simple elec- 
tron count with the usual valences for O 2- anions 
and La 3+ cations), a large crystal- field splitting in- 
side the Ni e g doublet is expected, with the d z 2 lying 
lower in energy. We can see this situation comparing 
the two band structures shown on the left of Fig. [2j 

It is not clear what spin-state one would expect in 
this compound, this would depend on the difference 
in energy between the intra-e^ crystal-field splitting 
and the Hund's rule coupling strength. These differ- 
ences can be derived from the two central panels of 
Fig. [2j Focusing only on the Ni atom (of which we 
have two in the structure) highlighted in Fig. [2] (the 
opposite spin would need to be used for the other 
one), if the down-spin d z 2 band were lower in en- 
ergy than the up-spin d x 2_ y 2, a low-spin state would 
be expected, and the d x 2_ y 2 up-spin band would be 
only half- filled leading to a metallic result. This low- 
spin state is the solution found in Ref. 20. However, 
our calculations yield a high-spin state to be more 
stable. In this case, as we see in Fig. [2j the up-spin 
d x 2_ y 2 is comparable in energy with the down-spin 
d z 2. However, there is a third energy to consider. 
We observe that the d z 2 bands come in bonding- 
antibonding pairs, due to the spatial confinement of 



the MO2 bilayers and the large coupling of those 
orbit als along the c-axis, similar to what happens 
in La4Nia08. That bonding- antibonding splitting 
leads to a gap opening around the Fermi level for 
this high-spin state that is favorable in energy ac- 
cording to our calculations. This result is consistent 
with the conductivity measurements performed in 
the compound. The fully occupied up-spin d x 2_ y 2 
stabilizes an in-plane AF ordering, according to our 
LDA+U calculations. 

A question remains as to why this system is not 
observed to have a long-range in-plane AF order- 
ing. We can try to compare the bandwidths of 
the d x 2_ y 2 band in the case of this compound and 
also in La4Nia08, for which we have calculations: 
the LasNi206 compound has a d x 2_ y 2 bandwidth 
of about 1.5 eV, whereas La4Nis08 has a band- 
width of only 1 eV. An increase in 50% in the band- 
width will lead to a significantly different behavior 
of the AF in-plane coupling. Also, in La2NiC>4, only 
short-range 33 AF order exists, which suggests that 
AF in-plane ordering caused by d x 2_ y 2 orbitals is a 
marginal occurrence in this series of compounds. 

To quantify the evolution of antiferromagnetism, 
the role of the O p orbitals should be considered. 
The fact that in cuprates AF is stronger— and sur- 
vives to larger values of doping could have to do 
with the O p bands playing a significant role in 
cuprates, but a lesser one in these layered nicke- 
lates. In the case of La3Ni2 06, the O p bands are 
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FIG. 2: Band structure of LasNi206, showing the fat bands for the different e^ states as labelled, calculated within 
the LDA+U scheme, with U= 4.75 eV for Ni. The d z 2 bands are bonding-antibonding split, while the d x 2_ y 2 bands 
show a larger bandwidth of about 1.5 eV. 



closer to the Fermi level (slightly more than 1 eV) 
than in I^NiaOg (around 2 eV below the Fermi en- 
ergy). Since the in-plane coupling is mediated by O 
p bands, differences of this kind will have implica- 
tions in the different magnetic properties observed 
in both compounds, where electronic structure cal- 
culations predict the same magnetic ground state, 
with a very similar molecular correlated insulating 
state produced by quantum confinement. 

This obvious difference (the lower position of the 
fully occupied O p bands compared to the cuprates) 
in the electronic structure between these molecular- 
insulating layered nickelates, both La3Ni2 06 and 
La4Ni308, 19 and the superconducting cuprates, 
where even Zhang-Rice singlets can occur— due to 
the position of the O p bands quite close to the Fermi 
level leads us to think that it could be interesting 
to explore the possibility of doping the oxygen sites 
with a larger isoelectronic anion like S. These S ions 
would create an anionic band with a larger band- 
width and situated higher in energy. The AF cou- 
pling in the plane could be enhanced as a result, 
and also the participation of the anion p bands in 
the electronic structure of these nickelates, making 
them more similar to the superconducting cuprates. 

Another difference lies in the fact that supercon- 
ducting cuprates have a very strong Jahn- Teller cou- 
pling that produces a large elongation of the oxygen 
octahedra. In the case of the insulating nickelates 
with an electronic structure closer to that of the 
Cu 2+ superconducting compounds, the environment 
is square planar. A much smaller coupling of the 
electronic degrees of freedom to the lattice would be 
expected in the nickelates as opposed to the case of 
the cuprates. 

One may ask: is it possible to have a pressure- 



induced metal-insulator transition in both I^NisOs 
and La3Ni206? The two energies controlling the 
spin-state in the Ni cations are the intra-e^ crys- 
tal field splitting, which can increase by reducing 
the a and b lattice parameters, and the bonding- 
antibonding splitting caused by the d 2 2-d z 2 strong 
bonding along the c-axis, which increases as the c 
parameter gets reduced. In a layered system like 
this, applying pressure will reduce the c parame- 
ter by a larger amount than the a, b in-plane lattice 
parameters. Pressure will move both d x 2_ y 2 bands 
higher in energy, and the antibonding part of the d z 2 
bands also higher in energy. It would be a question 
of how pressure affects the position of both bands 
that would determine a transition to a metallic low- 
spin state by applying pressure in this compound. 
With the experimental lattice parameters, we see 
from our calculations that the intra-e^ crystal field 
splitting between d z 2 and d x 2_ y 2 bands is about 2.8 
eV, whereas the bonding-antibonding splitting be- 
tween the d z 2 up-spin bands is only 0.8 eV. Thus, 
if applying pressure increases the crystal-field split- 
ting by a few tenths of an eV, a spin-state crossover 
may take place induced by pressure, leading to an 
energetically favored metallic state. 



V. LA3NI2O7 

The structure of LasNi207 25 is also layered, but 
compared to La3Ni2 06, all the apical oxygens neigh- 
boring Ni remain, leading to all Ni atoms being in an 
equivalent octahedral environment. A simple elec- 
tron count gives a Ni 2 5+ : d 7,5 state for both Ni 
cations. Experiments on the compound indicate a 
paramagnetic metal at all temperatures. However, 
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FIG. 3: Band structure of the La3Ni207 in the non-magnetic phase (left panel) and the spin- up and spin-down 
channels of the magnetic solution (central and right panel, respectively). The calculations were performed within the 
LDA+U scheme, with U= 4.75 eV in Ni atoms. In all cases, the Ni states are highlighted. The e g bands occupy 
states close to the Fermi energy, with dispersive d x i_ y i bands and relatively flat d z i close in energy. 



the formation of local moments and the appearance 
of dynamic AF coupling has been suggested by an- 
alyzing the susceptibility data. 35 Other studies con- 
firm the system to be a paramagnetic metal, show- 
ing signatures of a low-dimensional localized electron 
system* 8 - The existence of a charge density wave has 
been suggested in the past j 35 i 36 

Due to its complicated magnetic and electronic 
structure properties, we have studied two different 
solutions. One is a non-magnetic solution with an 
identical electronic structure in each spin channel 
(t2c/ 3 ec/ - 75 ), depicted in the left panel of Fig. [3j 
with the Ni states highlighted. Around the Fermi 
level, the e g states can be noticed, with the d x 2_ y 2 
band also crossing the Fermi level showing a 2 eV 
bandwidth. The d z 2 has a bonding and an anti- 
bonding branch due to the large interplanar coupling 
between those levels and the quantum confinement 
of the MO2 bilayers in the structure. This splitting 
leads to the d z 2 bonding states to be lower in energy 
and more occupied than the d x 2_ y 2 bands, which 
would otherwise be degenerate for the two Ni sub- 
layers. The bonding-antibonding splitting is slightly 
over 1 eV, larger than in La3Ni2 06, showing that 
a large hopping occurs also through an apical oxy- 
gen in the octahedral environment. The interplanar 
Ni-Ni distance is 3.1 A in La3Ni2 06 and 3.9 A in 
LasNi207, which might seem contradictory. It is 
however difficult to quantify the different effects of 
U (the on-site Coulomb repulsion introduced in our 
LDA+U calculations) here due to the different oc- 
cupation of that e g level compared to the case of 
La3Ni206- In both cases, we present results with 
the same value of U (4.75 eV). 



The other solution we have studied has a local 
magnetic moment for each Ni cation. This solution 
is lower in energy than the non-magnetic solution by 
a large 100 meV/Ni when LDA+U method is used 
with U= 4.75 eV. However, in such a metallic sys- 
tem the LDA+U description could be overestimating 
correlation effects, causing the energy difference to 
be too large. It is nevertheless interesting to study a 
magnetic solution, considering the experimental evi- 
dence for magnetic behavior. In such a solution, the 
Ni 2-5+ : d 7,5 cations have a spin value S=3/4, with 
the d z 2 orbital fully occupied, and the d x 2_ y 2 band 
half-filled for each spin channel. This quarter-filling 
of the in-plane e^ band would stabilize an in-plane 
AF coupling, if a localized electron picture is suffi- 
ciently valid in this itinerant compound. If correla- 
tion effects (quantified by U in our calculations) are 
eliminated from the calculations, the non-magnetic 
solution becomes more stable. The experimental ev- 
idence suggests correlated metallic behavior, which 
will require more sophisticated study. 



VI. LA3NI2O6.5 

The Mott metal-insulator transition in this series 
of compounds is expected to occur at the Ni 2+ com- 
pound. Starting to remove oxygen from the metallic 
La3Ni2 07, and using the usual valences for O and 
La, the transition should take place at La3Ni206.5, 
with Ni 2+ : d 8 cations, if these are in a high-spin S= 
1 state. The square-planar environment of the Ni 
cations in the insulating compound La3Ni206 sug- 
gests that on removing oxygen from La3Ni2 07, the 




FIG. 4: Band structure of LasNi206.5, showing the fat bands for the two different Ni cations in the structure. The 
square pyramidal Ni is on the left, the octahedral Ni is on the right. The majority spin figures show five occupied d 
bands. The minority spin figures show, above the Fermi energy, the flat d z i bonding-antibonding split bands below 
the d x 2_ v 2 band. 



first oxygens that can be removed are the apical oxy- 
gens in the NiC>6 octahedra. 

There are four inequivalent oxygens in the struc- 
ture of La3Ni2 07, the choice of an apical oxygen 
that leads to a ground state is more stable by sev- 
eral hundreds of meV/Ni compared to removing any 
other oxygen. Those total energies were calculated 
at U= 0, with a relaxation of the atomic coordi- 
nates done within the GGA scheme. The removal of 
an apical oxygen leads to two different environments 
for the Ni cations: a NiOe distorted octahedral envi- 
ronment and a square pyramidal MO5 environment. 
Once the choice of oxygens to be removed is clarified 
by total energy calculations, we have further relaxed 
the lowest energy state within the LDA+U scheme. 
Additional distortions of the oxygen cages surround- 
ing the Ni cations that occur lead to the opening of 
a gap at the Fermi level and the stabilization of a 
distinct S= 1 state in both the octahedral and the 
square pyramidal environments after relaxation. 

We can describe the slightly different electronic 
structures of the two inequivalent Ni atoms in the, 
respectively, square pyramidal and octahedral envi- 
ronments with the aid of the fat bands plots shown 
in Fig. [U Both cations are in an S=l state (t2 g 6 e g 2 )- 
The main difference between them is in the position 
of the d z 2 band. In the case of the octahedral en- 
vironment, it is in the same energy window as the 
d x 2_ y 2 band, particularly easy to see in the down 
spin channel on the third panel of Fig. [4j the flat 
band highlighted just below the Fermi level. In the 
case of the square pyramidal environment, the ab- 
sence of one of the apical oxygens leads the d z 2 band 
lying lower in energy than the d x 2_ y 2 band, due to 
the higher repulsion with oxygen undergone by the 



latter. The bandwidth of the band that transfers 
the in-plane AF coupling observed to be the ground 
state for this compound is small compared to that of 
the LasNi2 06, being only 0.6 eV. This small band- 
width could indicate that the in-plane AF ordering 
is even less likely in this borderline compound than 
in the more insulating ones, where it is also only 
marginally stable. Again, the O p bands that me- 
diate the coupling are several eV below the Fermi 
level, weakening the superexchange process. 



VII. SUMMARY 

In this paper we have presented a study of the evo- 
lution of the electronic structure in low- valence lay- 
ered nickelates with varying oxygen content (hence, 
Ni valence). The system La3Ni2 07_5, with 5= 0.0, 
0.5 and 1, has been studied. A metal-insulator tran- 
sition takes place between a molecular insulating 
state for 6= 1 to a paramagnetic metal and 5= 1, 
via a Mott-transition compound formed by Ni 2+ : d 8 
cations, with half- filled e g bands. The main differ- 
ence with layered cuprates lies in the absence of the 
O p bands close to the Fermi level; in nickelates, 
only Ni d states are close to the Fermi level. One 
consequence of this is a less strong in-plane AF cou- 
pling between Ni d x 2_ y 2 orbitals, and also a different 
range of existence of the AF ordering and the insu- 
lating phase with respect to oxygen doping. The 
insulating properties of LasNi2 06 are understood as 
a process due to quantum confinement of the MO2 
bilayer that leads to a bonding-antibonding splitting 
between the Ni d z 2 orbitals. The resulting gap gives 
the insulating behavior. 
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